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Consider any 1-parameter family of initial data such that data with parameter value p > p, form 
black holes, and data with p < p* do not. As p — ► p* from above ("critical collapse"), the black 
hole mass scales as M ~ (p — p») 7 , where the critical exponent 7 is the same for all such families 
of initial data. So far critical collapse has been investigated only for initial data with zero charge 
and zero angular momentum. Here we allow for 17(1) charge. In scalar electrodynamics coupled to 
gravity, with action R + \(d + iqA)4>\ 2 + F 2 , we consider initial data with spherical symmetry and 
nonvanishing charge. From dimensional analysis and a previous calculation of Lyapunov exponents, 
we predict that in critical collapse the black hole mass scales as M ~ (p — p*) 7 , and the black hole 
charge as Q ~ (p - p,) s , with 7 = 0.374 ± 0.001 (as for the real scalar field) and 5 = 0.883 ± 0.007. 
We conjecture that, where there is no mass gap, this behavior generalizes to other charged matter 
models, with 8 > 27. We suggest the existence of universality classes with respect to parameters 
such as q. 
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I. INTRODUCTION 

(N 

. Consider any 1-parameter family of initial data for general relativity such that data with parameter value p > p» 
' form black holes, and data with p < p* do not. Recently it has been discovered that as p — > p* from above ("critical 
collapse"), and unless there is a mass gap, the black hole mass scales as M ~ (p — p*) 7 , where the critical exponent 
\Q • 7 is the same for all such families of initial data. 

This behavior was first discovered for a real scalar field in spherical symmetry JlJ, has been observed in a variety 
of other matter models , and is strongly expected to occur in two larger classes |||7| , all in spherical symmetry. 
Critical phenomena were also found in axisymmetric pure gravity Q. Historically the second case to be discovered, 
it has so far remained alone in going beyond spherical symmetry, and in being vacuum. No counter-examples have 
O" 1 been found, except that when the field equations allow for a metastable static solution this solution creates a mass 
$_i , gap for certain initial data regimes [^) . 

The critical scaling of the mass can now be considered as understood |^||[l(| in terms of general relativity as a 
\ dynamical system, and of Lyapunov exponents, at least in spherical symmetry. For reasons of space we assume here 
that the reader is already familiar with these arguments. Our notation here is compatible with |j"o| . An introductory 
review is pd[| . 

Black holes are characterized by mass M, U(l) charge Q and angular momentum L, but so far critical collapse 
has been investigated only for initial data with zero charge and zero angular momentum. (Linear perturbations of 
extremal charged black holes (Q — M) were studied in fll2||.) Here, as a first step towards the generic situation, 
we allow for charge in the initial data and hence in the resulting black hole. We keep the restriction to spherical 
symmetry (and therefore exclude angular momentum). 
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II. SCALAR ELECTRODYNAMICS 



As a matter model with U{\) charge, we choose scalar electrodynamics coupled to gravity. This has the advantage 
of being a simple field theory, and of allowing arbitrary and independent initial distributions of energy-momentum 
and charge. Furthermore it is a generalization of two other models in which critical collapse has already been 
studied, namely the real and complex massless scalar fields (without electromagnetism). The action, in units where 
c = 2ttG = 1, is 

~i? - \g ab {V a + iqAa) 4> (V 6 - iqA b ) p ~ \g ah g cd F ac F bd , (1) 
where F a i, = V a A{, — VbA a . We restrict to spherical symmetry. For the metric we choose 
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ds 2 = -a(r, t) 2 dt 2 + a{r, t) 2 dr 2 + r 2 dVL 2 



(2) 



with remaining gauge condition a = 1 and regularity condition a — 1 (or a, r = 0) at the origin r = (before a 
singularity forms there). For the electromagnetic field we choose the gauge A r = 0, and At = at r = 0. 

In order to use dimensional analysis and scaling arguments later, we need the field equations in first-order and 
dimensionless form. For this purpose we introduce the dimensionless fields 



X = -£ ri Y = 
a 

With the notation Z ={X + ,X 



j> tt +iqA t (p), X±=X±Y, A = A t , F = —A t , r , g 

act 

, a, <?, A, F}, we can write the field equations formally as 
E(rZ, ri rZ, t , Z, qr) = 0, 



-At 



a 
a 



(3) 
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where the vector E of equations is polynomial in its arguments. (The matter stress-energy tensor and the matter 
equations of motions are manifestly polynomial because they have been derived from a polynomial matter Lagrangian. 
The components of the Ricci tensor can be brought into polynomial form after multiplication by certain metric 
coefficients.) Furthermore it is linear in the first two arguments, because the field equations are now first order. 

The arguments of E form a complete set of independent dimensionless quantities. To make dimensionless quantities 
out of q and the derivatives, we have combined them with r, taking advantage of the spherical symmetry. (The 
Ricci scalar R has dimension L~ 2 , and therefore F a \,F ah must have the same dimension. This means that A a is 
dimensionless. As the minimal coupling of matter to electromagnetism is via V a — > V a + iqA a , q has dimension L .) 

Before giving the equations explicitly, we introduce new dimensionless coordinates 



ln[l 



C 



&(r), 



(5) 



where r$ is an arbitrary fixed scale, and £o( T ) is a periodic function with period A. This choice of coordinates is adapted 
to the self-similar solution we shall be looking for Jl(|. In these coordinates, discrete self-similarity corresponds to all 
fields being periodic in r. (As a degenerate case, continuous self-similarity corresponds to all fields being independent 
ofr.) 

Rewriting the arguments of E in equation (Q), we obtain 



E(Z C , e c+ «°[Z, 



(l+QZ x ], Z, qr e T +t + t°) = 0. 



(6) 



By shifting the origin of r, one may set ro = q^ 1 to simplify the notation, but we want to stress here that q is always 
accompanied by e T in the field equations, and vice versa. This will be crucial later on. 
The complete field equations in these fields and coordinates, written out explicitly, are 



9X 

■Ac 

FX 



[1(1 



a 2 \X T \ 2 



2 F 2 ] X ± -X T ± [gX± tT + iqr e T (gAX ± + aF<j>)] 



= 1 -a 2 



2a 2 F 2 ) g 



= -a\l- a 2 
2 1 

= aX 



\XJ 



+ 2F 2 )] 



= -F 



+ iqr e T+c+ t a a- (cf>Y* - <j>*Y) 



A 



g^a 2 F 



r (C+&) l a a ff -i (| X+ |2 _ | X _ |2) + (1 + Q l fl {1 a 2 



IX, 



|A_ 



2F Z 



F, T = iqr^ag- 1 - (4>X* - 4>*X) + (1 + Q 



—F 



iqr e T+c+ t°a- (<pY* - <j)*Y) 



(l + QX + g^e-^+^Y 



iqroe 1 " A(j) 



(7) 

(8) 

(9) 
(10) 

(11) 
(12) 
(13) 

(14) 

(15) 



where £q = d^o/dr. The last three equations, which contain only r-derivatives, can be considered as constraints (at 
constant C) which are propagated in £ by the first six equations. 
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III. CRITICAL SOLUTION 



Critical collapse as we know it is dominated by a "critical solution" with two crucial properties: it is an attractor 
of co-dimension one, thus serving as an intermediate asymptotic, and it is self-similar, thus linking the large scale 
of the initial data with the small scale of the final black hole. The critical solution acts as an attractor within the 
hypersurface which separates black hole data from non-black hole data (the "critical surface"). It thus funnels all 
initial configurations on either side of the surface into two unique channels; one side goes to a black hole, and the 
other to dispersion. This explains the universality of the critical exponent. The power-law behavior of the black hole 
mass follows from the self-similarity of the critical solution. 

Finding a critical solution is a two-step process: one has to find a self-similar solution, and then check explicitly 
the number of its unstable perturbations. 



Self-similarity corresponds to periodicity of the solution in the coordinate r, while the field equations contain explicit 
factors of e T . Therefore no self-similar solution exists. The physical reason is that the coupling to electromagnetism 
introduces the scale q _1 into the action and field equations, thus excluding scale-invariant solutions. 

Does the absence of an exactly self-similar solution rule out critical phenomena? No, as we really only need a 
solution which approaches a self-similar spacetime as the latter approaches its singularity, that is in the limit of small 
spacetime scale. (Any spherical self-similar spacetime has a curvature singularity, with R ~ e~ 2r .) This leads us to 
the ansatz 



where each Z n (£, r) is periodic in r. As r — > — oo, Z is dominated by Zq, which is by assumption periodic in r. 
Therefore the form ( |l6| ) guarantees that the solution Z is asymptotically self-similar as r — > — oo. Like a self-similar 
solution, it has a singularity, at (t = 0,r = 0), which corresponds to r = — oo for any £. 

Now we argue that there are solutions of this form at least up to some distance from the singularity. We obtain 
equations for the Z n by substituting the ansatz into the field equations and separating powers of e T . Zq obeys an 
equation independent of the other Z n , while all the other Z n are coupled to Z„_i and lower. This allows us to 
determine the Z n recursively, starting with Zq. Therefore solutions Z of the full problem are mapped one-to-one to 
the solutions Zq of the zeroth order problem, which we shall see is simpler. 

We know from the study of self-similar solutions that Z is locally uniquely given by a nonlinear eigenvalue problem, 
arising from the boundary conditions of (i) periodicity in r and (ii) regularity at (a) r = (for t ^ 0, before the 
singularity forms) and (b) the past light cone of the singularity. By analogy, it is plausible that Z\ exists and is 
uniquely determined by the source term Zq and similar boundary conditions, and so on for all the Z n . 

Does the sum converge? The Z n are periodic in t. As £ — > — oo, which corresponds to r — > (for t ^ 0), the Z n are 
bounded by the regularity conditions we impose. £ — > oo corresponds to t — ► (for r ^ 0) and is only a coordinate 
singularity. The two examples for Zq discussed below have been continued, by a change of coordinates, up to the 
future light cone, and have been verified explicitly to be bounded. As the higher Z n obey equations similar in type 
to Zq, we assume that all Z„((, t) are bounded. As for the growth of the Z n with n, it seems very likely that it is 
bounded as Z n < Ae~ Tin for some n, so that the solution converges for r < t\. The fact that it does not converge 
for all r is not crucial to its usefulness. Even when the critical solution is exactly self-similar and therefore known to 
exist for all r, an actual collapse spacetime is not even approximately self-similar outside a bounded spacetime region. 

What is the meaning of the expansion (|l^) in powers of e r ? As e T always appears together with q in the field 
equations, it is also an expansion in powers of g, that is in the coupling of the matter to electromagnetism. For example, 
while Aq vanishes, A\ is the electromagnetic field created by 4>q, and 4>2 is the reaction of 4> to that electromagnetic 
field, and so on. In particular, the equation for Zq is obtained by writing Zq for Z in the field equations while setting 
all appearances of q equal to zero. In effect, Zq obeys the field equations for a complex scalar field coupled only to 
gravity, but not to electromagnetism, which, as we anticipated, is a simpler problem. Exactly two self-similar solutions 
of that problem, or candidates for Zq, are known: 

1) The critical solution found by Choptuik [p] JT^ , ^0| ] has no charge, because the scalar field is real (up to the 
remaining U(l) gauge freedom <f> — > e lc <fi with c a constant). This means that all the Z n for n > 1 have vanishing 
source terms and therefore vanish, or in other words that the Choptuik solution is already a solution of the full 
problem. (Strictly speaking, we have to complete the real Choptuik solution by a vanishing A and F and a real <p 



A. Asymptotically self-similar solutions 



oo 




(16) 
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calculated from (15). <fi is periodic because in the Choptuik solution X and Y have vanishing mean value. This is 
essential, but was not automatically guaranteed.) 

2) The continuously self-similar, charged solution found by Hirschmann and Eardley 0] does have charge. This 
means that it gives rise to a solution Z of the full problem, including electromagnetism, only when all the terms Z n 
are added up. 



B. Perturbation spectrum 



Now we come to the second criterion for a solution to be a critical solution, namely the presence of exactly one 
linear perturbation mode that is growing as r — > — oo. Anticipating the next section, we note that the general linear 
perturbations SiZ of Z (i labels the independent modes) have a similar series form to the background critical solution 
Z itself. As for the background expansion(whcre Z is completely determined by Zq), 5iZ is specified by giving its 
leading term S^Z . This obeys an equation where (a) once more the terms with q have been set equal to zero, and 
where (b) the background dependence is only through Zq. This means 5toZ on its own is a complete perturbation, in 
the theory without electromagnetism, of Zq, which is itself a complete background solution of that theory. 

The spectrum of the perturbations is determined, in a linear eigenvalue problem, together with the S^Z (and is 
not influenced by the higher 5i n Z). In consequence, the spectrum is the same in the presence of the coupling to 
electromagnetism, or in its absence. We can therefore recycle earlier work, where the perturbation spectra of the 
Choptuik and Hirschmann and Eardley solutions were determined in the theory without electromagnetism: their 
spectrum in scalar electrodynamics is just the same. As in particular the number of unstable modes is the same, we 
immediately conclude that the Choptuik solution is an attractor of co-dimension one @, that is a genuine critical 
solution, even in the full theory with electromagnetism, while the complex solution based on (but not identical with!) 
the Hirschmann and Eardley solution has three unstable modes (To state it once more, its perturbation mode 

functions are different in the full theory, but their Lyapunov exponents are the same.) 



C. Uniqueness of the critical solution 

Collapse simulations Jl6| suggest that the Choptuik solution is in fact the only critical solution, and a global 
intermediate attractor, for the free complex scalar field not coupled to electromagnetism. This type of evidence can 
never be complete, because the entire phase space can never be probed. It only suggests that any other other critical 
solutions have rather small basins of attraction. As we have seen, self-similar solutions and their perturbations are 
mapped one-to-one between the free complex scalar field and the theory with electromagnetism. This means that 
critical solutions are mapped one-to-one between the theories. Their basins of attractions, however, could be very 
different in size, so that the counterpart of a hypothetical critical solution for the free complex scalar, with a basin of 
attraction so small that it would have been overlooked so far, could play an important role in scalar electrodynamics. 

We must also consider the possibility of critical solutions not of the form (|l6|). Two such possibilities have occurred 
to us, but can be definitely ruled out. 

Any critical solution of the form ( |l6| ) describes a situation where charge becomes less and less important in the 
final black hole as one fine-tunes along some one-parameter family of initial data with the aim of making black holes 
of ever smaller mass. (This does not exclude initial data with an appreciable charge-to-mass radius, but most of 
that charge must be radiated away along with most of the mass.) The other possibility for the charge allowed by 
the limit Q < M for black holes is that Q/M approaches a constant as M — > in the process of fine-tuning. This 
is excluded: The black hole charge Q for a given scalar field evolution must change sign when q changes sign. As q 
always appears in the company of e r , Q must be suppressed with respect to M by at least one power of e T as small 
scales are approached. 

A second possibility (suggested by ||) is that of another critical solution which is not self-similar but static, 
giving rise to a mass gap in the one-parameter families of data. (We come back to this possibility below, when we 
consider other matter models with charge.) But massless scalar electrodynamics has no static solutions, even unstable. 



"Charged boson stars" exist only if one adds a mass term m 2 </> 2 to the action, with AnGm 2 > q 2 |17 
We come to the surprising conclusion that the Choptuik solution is a critical solution in scalar electrodynamics, 

and, as far as we can see, the only one, even though it has no charge itself. (In particular, the Hirschmann and 

Eardley solution, although charged, is not a critical solution.) 

On goin g work of C. G. shows that the critical solution for spherical SU(2) Einstein- Yang-Mills collapse [9 is of the 

form (|lq), with all the terms Z n nonvanishing. This indicates that if they vanish for n > in scalar electrodynamics, 

this is accidental. Details will be published elsewhere. 
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IV. LINEAR PERTURBATIONS 



We now consider the real critical solution of Choptuik and its linear perturbations. Because the background is real, 
these split naturally into two kinds. The first are perturbations 5Z with S(j) purely real. In the following we call these 
the "real" perturbations. The general perturbation of this kind is of the form 

oo 

SZ ieal (C,r) =J2°i e XiT SiZ(C,r), (17) 
»=i 

where each SiZ is periodic in r with period A, and where the Ci are free parameters. They were already calculated 
within the real scalar field model, in In particular there is precisely one eigenvalue A with 3?A < in the spectrum 
{Xi}, namely Ai ~ —2.674. 

The linearized equations for the remaining, not purely real, perturbations are of the form 

SZ x =ASZ + BSZ^ T ,+e T CSZ, (18) 

where the dimensionless coefficients A, B and C are periodic in r. The "real" perturbations obey the same equations, 
but the term proportional to e T (or, as always, proportional to q) vanishes for the real perturbations because they 
don't carry charge. Its presence means that the 5iZ cannot be periodic in r. Instead we have 

oo 

^ imaginary (C r) = J2 A e^ T 5iZ(C, r), (19) 

z=l 

where each independent mode SiZ is now expanded as 

oo 

S t Z(C,r) = Y / e nT S m Z(C,r), (20) 

71 = 

where only the 5i n Z are periodic in r. This expansion is exactly analogous to that of the asymptotically self-similar 



background solution, with the difference that, due to the linearity, the ansatz ( |19| , |20| ) is generic, and (16) is not. The 
5i n Z obey the coupled equations 

S M Z X = (A + fHB) S i0 Z + B 6 i0 Z, T , (21) 
S m Z x = (A + ^B + nB) S m Z + B 5 in Z, T + C 5 in - X Z, n > 1. (22) 

The equation for S^Z describes purely imaginary perturbations of the scalar field, not coupled to the electromagnetic 
field or the metric. It is homogeneous in S^Z ', and is complemented by periodic boundary conditions (with period A) 
in t, and by regularity conditions at Q = — oo and £ = 0. A solution exists only for discrete values of fi. From this 
equation one obtains the spectrum {fa} of perturbations. This equation is in fact that for the perturbations, around 
the real solution, of the complex scalar field model without electromagnetism. They have already been considered in 
[fTo| , with the result that $lfa > for all i. 

The terms Si n Z for n > 1 are determined both by boundary conditions and by the source terms 8i n ^\Z . As the value 
of fa has already been fixed as an eigenvalue in the equation for S^Z , the boundary conditions admit no homogeneous 
solution for n > 1, and the solution for n > 1 is proportional to the source term bi n ^\Z. The electromagnetic field 
comes in to order e r , that is q, because the purely imaginary perturbation 5io<p of the real background scalar field 
</>* gives rise to a charge distribution and therefore generates an electromagnetic field. This charge distribution is 
ultimately responsible for the charge of the final black hole. To order e 2r , or q 2 , the scalar field perturbations acquire 
a real part and in consequence couple to metric perturbations, thus also creating a gravitational field. Nevertheless, 
we shall refer to these perturbations as the "imaginary" perturbations, because to leading order they are purely 
imaginary perturbations of the scalar field. 

V. MASS SCALING 

For the generic solution close to the real self-similar solution Z*, we now have the form 

oo oo oo 

Z(C,r)~Z»(C,T)+^C 4 (p) e x > T *Z(C,t)+£A(p) ]T e" T S m Z((,r), (23) 

i—1 i—1 n—0 



5 



where the first term is the critical solution, and the second and third terms its "real" and "imaginary" perturbations, 
as discussed above. The amplitudes and Di of the perturbations depend on the initial data in general and hence 
on the parameter p of a given one-parameter family of initial data in particular. 

As r — * — oo, we can neglect all perturbations but the one growing mode, associated with Ai. As we are interested in 
the asymptotic behavior of the charge however, we also keep the most slowly decaying of the imaginary perturbations 
(which alone carry charge) , associated with fii . By definition we obtain the precisely critical solution for p = , and 
so we must have Ci(p*) = 0. Expanding Ci(p) and D\(p) to leading order, we obtain 



Z(C, t) ~ Z*(C, r) + (p~ p.) e AlT SiZiC, t) + D x (p.) £ e nT 6 ln Z((, r). (24) 

To keep the notation compact, we define, following ]Io| , 

r . F>r„ 1 

to = n + (25) 



r*(p) = In ( - — — ] , n = -i In 
Ai \ p* J Ai 



-i 9d : 
amp 



where e is an arbitrary small constant (independent of p) . If we now fix r = tq in the approximate solution ( p4[ ) , we 
obtain a p-dependent family of Cauchy data, namely 

Z p (r) = Z{(, t ) = Z* (ln^, T ^j + e S,Z fin ^, r ^) + K(p) £ e nT ° 5 ln Z (ln^, t q ) (26) 

where 



n=0 



p = r e T "+^ T "\ K{p) = Difay™ . (27) 



Kip) is small, even compared to e, if e r ° is sufficiently small, that is, if (p — p*) is sufficiently small. We therefore 
treat the terms proportional to K(p) as a linear perturbation throughout. 

But now we consider the exact, nonlinear evolution of the data Z* + t8\Z , without treating eS\Z as a perturbation 
any longer. This is necessary because its presence makes a qualitative difference at late times. If it has one sign, a 
black hole is formed. If it has the other, the matter disperses. Because the solution at late times is no longer even 
approximately self-similar, we go back to the coordinates r and t. 

The data Z* + eSiZ are purely real, and evolve to a purely real solution, with vanishing electromagnetic field. 
Consequently, the equations determining the solution do not contain the term qr, and are scale-invariant. Therefore 
the entire solution depends on r p in the simple way 

r t — t 

Z(r,t) = f(r,t,T ), where f=— and t = E , (28) 

Tp Tp 

where the (irrelevant) shift in t is t p = r^e TQ . f obeys the equation 

E {fir, fit, /, 0) = 0, (29) 

with initial data 

/(f, i=0,T ) = Z* (lnf, r ) + e SiZ (lnf, r ) . (30) 

In particular we know that the mass of the black hole that forms in this solution must be a multiple of the underlying 
scale r p , namely 



.M(ro) ; ( 31 ) 



where /i(r) is a function that we do not know, but which is periodic with period A in r. The black hole mass as a 
function of the family of initial data and the parameter value p is then 



M(p) = r l^—^-j e/i+»lMf>)+ri) ; (32) 

where v{t) = /u(r) + £o(t) is a universal wiggle superimposed on the basic power-law behavior, and the constant n 
(defined above) depends on the initial data family. 
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VI. CHARGE SCALING 



So far we have only repeated |T^] (which itself is a generalization of Now we consider the effect of the 

perturbation proportional to K(p) in the initial data for /. This linear perturbation, say Sf, obeys the linearization 
of equation (^), which is of the form 

r Sf, t = Ar Sf, r + B Sf + qrC Sf. (33) 
Rewriting this in terms of f and f, we obtain (using the definition ( ]27j ) of r p ) 

f 8f- t = Af 6f, f + B 5f + qr e To+i " {To) Cf Sf, (34) 

where the dimensionless coefficients A, B and C are functions only of f and i (as well as of the parameter tq 
characterizing the background solution /). To obtain a solution for all small values of the parameter e T ° at once, we 
expand in powers of it. For the critical solution plus linear perturbation we then have 

oo 

Z P (r, t) = f (f , t, T ) + K{p) e nT ° S n f (f, i, r ) . (35) 

The expansion coefficients S n f obey the equations 

f S f t t = Af 6of, f + B Sof, (36) 
f Snf t t = Af 6 n f t r + B5 n f + qr e£°™ Cf 5 n ^f, n > 1, (37) 

with initial data given by 

6 n f(f, i=0,T O )= S ln Z (In f , r ) . (38) 

The perturbation Sf gives rise to a perturbation of the black hole mass M, but we ignore this here as a subdominant 
effect. We are however interested in the black hole charge Q, which only comes in through Sf. We do not need to 
calculate Sf to see how Q scales. It is sufficient to note that the charge-to-mass ratio is dimensionless and must be 
odd in the coupling constant q, and hence e T °, and therefore must go as 



1 (To)+e3To (S) 3 (To)+o ( e5ro ) 



(39) 



Taking only the dominant term, and putting it all together, we obtain 



Q(p) = M(p) K(p) e T ° [%) (t ). (40) 



M 
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After regrouping terms, we obtain the final result 



Q(p)=r (^-^J 1 e T2+ ^ Tl+T '^\ (41) 

where ti = (pi + 2)ti + lnl?i(p») is a new family-dependent constant, t\ is the same family-dependent constant as 
before, and 7t(t) = v(j) + In (Q/M)i(t) is a new universal wiggle. 

Numerical values for the Lyapunov exponents are Ai = —2.674 ± 0.009 [|l0| and p\ — 0.362 ± 0.012. p\ has been 
calculated by the same method [ p~0| as Ai, but converges somewhat more slowly with decreasing step size, so that 
the estimated numerical error in fi\, and in consequence in 5, is somewhat larger. We obtain critical exponents 
7 = -1/Ai = 0.374 ± 0.001 for the mass and 5 = -(m + 2)/Ai = 0.883 ± 0.007 for the charge. 



VII. OTHER CHARGED MATTER MODELS 



Given that our arguments did not rely on the exact form of the field equations, we can generalize them to any 
matter model in spherical symmetry where the matter is coupled to electromagnetism only via the J7(l)-covariant 
derivative D a = V + iqA a . Then it follows from dimensional analysis that whenever one casts the field equations in 
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first-order, dimensionless form, they must be of the form (|f|). Z would stand for another set of fields, and E need not 
be a polynomial. The argument showing that the electromagnetic interaction can be neglected asymptotically in the 
strong field/ small scale regime of the critical solution would go through as before, even if the critical solution itself 
carries charge. (Strictly speaking, we require that a polynomial form of the field equations exists for our arguments 
to apply directly, but we hope that this technical requirement can be relaxed.) Now we need to consider three cases 
separately: 

(1) The critical solution does not carry charge, only its perturbations do, which brings in Lyapunov exponents. Let 
Ai be the one eigenvalue with -RA < 0. 

(la) The unstable mode itself carries no charge. This case is similar to scalar electrodynamics. Let fix be the 
eigenvalue associated with charge which has the smallest real part. We must have > 0, because the critical 
solution, by definition, has only one unstable mode. From 7 = — 1/Aj and 5 = — (/J-i +2)/Ai we then find the relation 



(lb) The unstable mode itself carries charge. In this case, we obtain the charge, as well as the mass, from the 
nonlinear evolution of the data + eSiZ, neglecting all the other perturbations. The solution arising from these data 
is no longer of the simple form (p8|), due to the presence of qr in the equations (^). But we can obtain the solution 
for all small r p by expanding in qr p , using the fact that qr = fqr p . The solution (^8|) is now replaced by the more 
general 



The mass M is again proportional to r p , but this time only to leading order in qr p (or e T °). The charge-to- mass ratio 
is now Q/M = (qr p )(Q / 'M)\{tq) + 0(qr p ) 3 , without the prefactor K(p), so that the charge to leading order is simply 
proportional to r p , or S — 27, with 7 = — 1/Aj.. 

(2) Finally, the critical solution itself may carry charge. Then the higher terms in the expansion ( |l6| ) do not vanish 
identically. This would have been the case for example for the solution fL4| . There the charge would have appeared 
to order e T of the background expansion (|l6|). We should still need to consider the data Zq + eSiZ, where 8\Z is the 
one growing perturbation (charged or not), in order to determine r p , the spacetime scale on which the solution leaves 
the intermediate attractor. But the electromagnetic field giving rise to the black hole charge Q is now dominated by 
Ai (the component A of Z\\ which is once more down a factor of r p from Z , so that Q is again proportional to r^, 
or S = 27. 

We must make one proviso, namely that a mass gap may exist in critical collapse. Recent work on critical collapse 
of the Einstein- Yang-Mills system has illustrated this [||: If the model admits a static or oscillating solution with 
precisely one unstable mode (such as the Bartnik-McKinnon solution in Einstein- Yang-Mills) , this solution acts as an 
alternative intermediate attractor. As this intermediate asymptotic is not self-similar, but instead has a finite mass, 
it gives rise to a mass gap in some region of initial data space. Our results then only hold for that part of initial 
data space where formation of the static solution is avoided. For a mass gap to occur when a critical solution also 
exists, it is essential that a static or oscillating solution not only exists, but has exactly one unstable mode. If it 
was unconditionally stable, there would be a third possible outcome besides black hole formation and dispersion - 
formation of a star. (This is of course the situation in astrophysics.) If it had more than one unstable mode, it would 
be a lesser attractor than the critical solution and would be "missed" by almost all one-parameter families of initial 
data. 



Finally, we note a consequence of our work that does not concern charge. By exactly the same argument which 
shows that in the critical solution with charge electromagnetism can be neglected asymptotically, one can show that 
the mass term can be neglected asymptotically for the massive scalar field. (One only has to replace qr by mr in 
equation (|j) and the following equations derived from it.) As the critical exponent for the black hole mass is unchanged 
by the coupling to electromagnetism, so it is, by the same argument, when a mass term or a more general polynomial 
scalar field self-interaction is added to the action. (This was already known flq| .) 

Our argument generalizes this to any parameter m of dimension (length)" 1 in geometrical units G = c = 1 appearing 
in a polynomial form of the field equations. (Note that as we insist on the equations being polynomial in the fields Z 
and parameter m, we implicitly also allow parameters of dimension (length) - ™, but only for positive, integer n.) 

In the language of critical phenomena, theories with different values of such a dimcnsionful parameter form "uni- 
versality classes" pT| . The analogy with critical phenomena in statistical mechanics seems good enough to use this 
term deliberately. 



6 > 2 7 . 



OO 




(42) 



VIII. UNIVERSALITY CLASSES AND THE RENORMALISATION GROUP 
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The renormalisation group in statistical mechanics acts on the phase space by blocking degrees of freedom, and a 
change of scale. If one demands that the partition function remain invariant, this induces an equivalent action on (the 
parameters of) the Hamiltonian p9| . The equivalents of the spins and their Hamiltonian in statistical mechanics are 
the fields Z(r) and their equations of motion in critical collapse. The renormalisation group in critical collapse acts 
on the phase space by a time evolution (in t) followed by change of scale (in r ), as Z(r,t) -► Z{e- A r,e- A t). (One 
can combine the time evolution and rescaling into a time evolution in r, at constant £, using the freedom of lapse and 
shift in general relativity. The critical solution is by definition invariant under this transformation.) If one demands 
that the rescaled data evolve in the same way, this action induces an equivalent action on (the parameters of) the 
field equations. 

In the case of scalar electrodynamics, this action is rather simple: In the field equations, q and r appear together 
as qr. If one changes the scale as r — > e~ A r, this has the same effect as q — > e~ A q. 

Critical phenomena in gravitational collapse correspond to very small scales (compared to the scale of the initial 
data and the scale q^ 1 ), in contrast to critical phenomena in statistical mechanics, where they correspond to very 
large scales (compared to the scale of the microscopic physics). In the limit of small scales, the parameter q becomes 
"irrelevant" , to borrow another term from statistical mechanics. 

If there are two (or more) parameters such as to or q, universality does not hold, as mxjmi is a dimensionless 
parameter, which generically has qualitative effects. In massive scalar electrodynamics, for example, the value of 
| q | /to determines if static solutions ("charged Boson stars") exist. Universality may be recovered if the two scales are 
very different, that is for m\/rri2 — > or TO2/TO4 — > 0, but this limit need not be regular. 



IX. CONCLUSIONS 



We predict that in critical collapse of massless scalar electrodynamics in spherical symmetry, the black hole mass 
scales as M ~ (p — p*) 7 , and the black hole charge as Q ~ (p — p*) s , each overlaid with a universal wiggle, with 
7 = 0.374 ± 0.001 (as for the real scalar field) and 5 = 0.883 ± 0.007. 

We have gone beyond the restriction to uncharged initial data (and hence black holes), but have kept the restriction 
to vanishing angular momentum. We predict the appearance of a new critical exponent 5, and give its numerical 
value. Verification of our predictions in numerical collapse simulations should be straightforward. 

Furthermore, we predict that in spherical critical collapse of matter models with minimal coupling V — > V + iqA 
to electromagnetism, both mass and charge scale as universal power laws (times a universal wiggle), with 6 > 2j, 
so that the black hole charge always disappears faster than the black hole mass. Depending on the matter (not in 
massless scalar electrodynamics) this behavior may hold only for parts of the initial data space, with a mass gap in 
other parts. 

We also find that adding terms with a parameter to of dimension (length) _1 in units c = G = 1 (such as a scalar 
field mass, or the electromagnetic coupling in the present paper) - in a form of the field equations which is polynomial 
in the fields Z and parameter to - does not change the critical exponent for the black hole mass. In this sense models 
with different values of such a parameter form universality classes. 

Notes added: After this paper had been submitted, our value of the critical exponent for the charge was confirmed 
in numerical collapse situations by Hod and Piran |20|. Universality classes in critical collapse were independently 
described by Hara, Koike and Adachi |2l|. 
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